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The one-dimensional extended bosonic Hubbard model has been shown to exhibit a variety of 
phases ranging from Mott insulator and superfluid to exotic supersolids and Haldane insulators 
depending on the filling and the relative value of the contact (U) and near neighbor (V) interaction 
strengths. In this paper we use the density matrix renormalization group and the time evolving 
block decimation numerical methods to study in detail the dynamics and excitation spectra of this 
model in its various phases. In particular, we study in detail the behavior of the charge and neutral 
gaps which characterize the Mott, charge density and Haldane insulating phases. We also show that 
in addition to the gapless modes at k — 0, the supersolid phase exhibits gapless modes at a finite k 
which depends on the filling. 


I. INTRODUCTION 

The bosonic Hubbard model (BHM) has continued 
to attract interest since its introduction by Fisher et 
al. pQ. This interest stems from its use to understand 
many physical phenomena such as the effect of disor¬ 
der on superfluids and the appearance of the compress¬ 
ible Bose glass phase [I], quantum phase transitions be¬ 
tween strongly correlated exotic phases etc. Interest in 
the BHM intensified with the experimental realization of 
Bose-Einstein condensates and the ability to load them 
in optical lattices |2j. Under experimentally realizable 
conditions, these systems are described by the BHM and 
its extensions [3] with highly tunable parameters and in 
one, two and three dimensions which makes them ideal 
for studying quantum phase transitions and exotic phases 
in strongly correlated systems. 

Increasingly, over the last several years, the physics of 
strongly correlated quantum systems has focused on the 
existence and properties of unconventional phases and 
phase transitions. In addition to well studied Mott in¬ 
sulating behavior caused by strong on-site repulsion at 
commensurate filling, extensive quantum Monte Carlo 
(QMC) simulations have shown that a strong enough 
near neighbor repulsion can lead to insulating incom¬ 
pressible density wave order (CDW) at integer and half 
odd integer fillings. Topological phases, such as the Hal¬ 
dane insulator which is a gapped phase characterized by 
a non-local (string) order parameter [4, 5] can be also 
found in the extended BHM [SHE] in one dimension. Fi¬ 
nally, doping these phases can lead to phase separation 
or to supersolid (SS) phases [9H22] . 

Even though the phase diagram of the extended BHM 
is now well understood, the excitation spectra of the var¬ 
ious ground states have been less studied [23l [24] , es¬ 
sentially because the numerical methods providing the 
ground state properties, such as exact diagonalization or 


QMC, become limited in the dynamical domain. More re¬ 
cently, for quasi-ID systems, the extension of the density 
matrix renormalization group method (DMRG) to the 
time domain or, equivalently, the time evolving density 
matrix method (TEBD) have proved to be extremely suc¬ 
cessful in probing the dynamical properties of the system, 
thereby providing reliable excitation spectrum [[25l - f27] . 
In this paper, we extend our work in [Si to study the ex¬ 
citation spectrum of the one dimensional extended BHM 
in different phases, namely the Mott insulating phase 
(MI), the Haldane insulating phase (HI), the charge den¬ 
sity wave phase (CDW) and the supersolid phase (SS). 


The paper is organized as follows. In section [TT] we 
present the model and the various methods to compute 
the ground state properties and excitation spectrum. In 
section m we present the dependence of the dynamical 
structure factor on the near neighbor repulsion, V, at 
fixed filling, n m 1, and fixed value of the contact re¬ 
pulsion, U and hopping, t. In section [IV| we study the 
dynamical structure factor in the SS phase for different 
fillings and explain its main properties using a mapping 
of the extended BHM to the Heisenberg model for a spin 
1/2 chain in a finite magnetic field. In addition, we de¬ 
scribe the evolution of the spectrum across the SS-SF 
transition. In section [V| we discuss the qualitative dif¬ 
ferences found in the dynamical structure factor in the 
phase which is obtained by underdoping the half-filled 
CDW and in the SS phase. A summary of results and 
conclusions is in section IVll 
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II. MODEL AND METHODS 
A. The model Hamiltonian 

The one dimensional extended BHM we shall study is 
described by the Hamiltonian, 

H = -t^2(a\a i+ 1 + aj +iai ) + ^ (ra* - 1) 

i i 

+V'Y^n i n i+ i. (1) 

The sum over i extends over the L sites of the lattice, 
periodic boundary conditions were used in the QMC sim¬ 
ulations [28] and open conditions in the DMRG and the 
TEBD simulations. The onsite repulsive interaction en¬ 
ergy, [/, is put equal to unity and sets the energy scale. 
The time scale being given by h/U , it is also put equal to 
unity. The operator a i (aj) destroys (creates) a boson on 
site i, rti = a\a i is the number operator on site i, t is the 
hopping amplitude and V is the near neighbor repulsive 
interaction parameter. Since we will typically study the 
system in the canonical ensemble, we did not include a 
chemical potential term in H . 

The charge gap is given by, 

A c (n) = fi(n) - fjb(n - 1) (2) 

= E 0 {n + 1) + E 0 {n - 1) - 2 E 0 (n) (3) 

where the chemical potential is given by /i(n) = Eo(n + 
l)—Eo(n) and Eq{u) is the ground state energy of the sys¬ 
tem with n particles and is obtained both with QMC and 
DMRG. The neutral gap, A n , is obtained using DMRG 
by targeting the lowest excitation with the same num¬ 
ber of bosons. For the CDW and HI phases, the chem¬ 
ical potentials at both ends are set to (opposite) large 
enough values, when using DMRG, such that the ground 
state degeneracy and the low energy edge excitations are 
lifted [6[ 0 [24j |29] . 

For a bosonic filling n close to unity, the Bose-Hubbard 
model can be reasonably approximated by the AF spin-1 
Heisenberg model: 

H S = E J(S?S* +1 +S?S? +1 )+\S?S? +1 +D (S?) 2 , (4) 

i 

where A is the axial anisotropy and D the ion anisotropy. 
One has the following mapping between the parame¬ 
ters [7]: J = —tn , A = V and D — U j 2. 

B. Time Evolving Block Decimation 

As mentioned above, the excitation spectra are ob¬ 
tained using the TEBD, first in imaginary time to ob¬ 
tain the ground state, then in real time to compute 
the density-density correlation function. In each case, 
we have used a number preserving algorithm. We have 


checked that the ground state properties (energy, site 
density, double occupancy...) obtained with the TEBD 
exactly match the properties of the ground state obtained 
from the DMRG, using the ALPS library m- The space 
and time correlation functions ( Ai(T)Bj ), where (• • •) is 
the ground state average and where Ai(T) is the time 
evolution of the operator Ai in the Heisenberg picture, 
have been obtained by writing, 

(Ai(T)Bj) = e iEGsT (GS\Aie~ iHT Bj\GS), (5) 

where Eqs is the ground state energy. Therefore, com¬ 
puting the correlation function can be done as follows: 

• From the matrix product state (MPS) representa¬ 
tion of the ground state, | GS), one computes the 
MPS of the initial state |\£(0)) = Bj\GS). 

• The state |\I/(0)) is evolved using the real time 
TEBD, providing the MPS of |^(T)), and thereby 
allowing the computation of the correlation func¬ 
tion as {GS\Ai\V(T)). 

In what follows, we focus on the density-density cor¬ 
relation, i.e. Ai = Bi = rii — (n*), more precisely, the 
initial state consists of creating a density excitation in 
the middle of the chain: Bj\GS) with j = L/ 2. We 
then compute (Ai(T)Bj) for all sites and times T up to 
150, with a time step equal to 0.1 (a smaller timestep 
was actually used for the propagation). Finally, the dy¬ 
namical structure factor S(k,uj) is computed from the 
Fourier transform of the density-density correlation with 
respect to i—j and T. In order to smooth out oscillations 
caused by the finite time window, we actually compute 
the Fourier transform of (Ai(T)Bj) exp (—4T 2 /Z^ ax ). 


III. MOTT-HALDANE-CDW TRANSITIONS 

The qualitative description of the different phases is 
based on the Heisenberg model EQ. ©• More precisely, 
the phases are characterized by the values of the string 
order parameters: 

O a = lim (S“e" E S« s p5“). (6) 

\i-j\^oo 

characterizing a loose antiferromagnetic order along the 
different axes a = x,y,z. They are associated with an 
underlying non-local discrete Z 2 G Z 2 symmetry of the 
Heisenberg model, Eq.Q [30]. In the large-D phase (the 
MI state for bosons), the O a vanish. In the Haldane 
phase, the two discrete Z 2 symmetries are broken, re¬ 
sulting in nonvanishing string order O a . Finally, in the 
Ising phase (the CDW for bosons), only the Z 2 symmetry 
along the z-axis is broken such that only the string order 
O z is non-vanishing. Note that in that phase, the string 
order and the antiferromagnetic order (the CDW order 
for bosons) are equivalent [34] . 
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The lowest elementary excitations from the ground 
state (ka = 0, S z = 0) are part of a triplet, one neu¬ 
tral e(°\k) (SN = S z = 0,), two charge ones e (±) (k) 
(SN = S z = ±1) [23[ |35], where SN corresponds to 
change in the total number of bosons. In each sector, 
one defines a gap which corresponds to the minimum of 
the elementary excitations over all k values: and 


The minima are located either at ka = 0 or ka = tt. 
From linear response theory, the structure factor reads 

| (fc, m\Sh\GS )\ 2 


S(fc, uS) oc 


w + ir) + - E\ 


k,m 


( 7 ) 


where |fc,m) denotes the different excited states of H 
for a given momentum /c, and E^m the corresponding 
energy. For single excitations, one simply has = 

e^ rn \k). Doubly excited states for fixed k are made of 
two single excitations: \q, m\ k — g, m'), corresponding to 
an energy e ( ^ rn \q) + \k — q). 

By definition, the charge gap of the system is A c = 
+ G^~\ i.e. the minimum energy for adding a par¬ 
ticle plus the minimum energy for removing a particle 
(increasing or decreasing S z , in the Heisenberg model). 
The neutral gap corresponds to the minimum of ei¬ 
ther the elementary neutral excitations, i.e. G^\ or of 
e^ + \k — q) + e^~\q), i.e., a combination of two charge 
excitations. Since, the minimum of both the elementary 
charge excitation e^(q) is attained at either q = 0 or 
q = 7r, the minimum of the two-particle excitation neces¬ 
sarily takes place at k = 0 and has the value G^ -\-G^~\ 
and corresponds then to the lower bound of a two-particle 
continuum. In short, the neutral gap value is given by 
the minimum of G^ and G^ + G^~\ 



Figure 1. Mott - Haldane - CDW transitions at fixed t/U = 
0.25. Around Vc ~ 0.75F, the neutral and charge gaps start 
to differ, indicating that G^O) < + G G); the gap for the 

single particle neutral excitations is smaller than the gap for 
the two particle excitations. The fact that only the neutral 
gap vanishes at the HI-CDW transition is a signature that 
the string order O z along the z -axis remains finite, where O x 
and O y orders vanish. 

It is well-known that in the Haldane phase, the neutral 


gap changes from one type to the other mmm and 
is emphasized in Fig. [lj where around Vc ~ 0.751/, the 
neutral and the charge gaps start having different values. 
For lower V values, one has A n = A c = G^ + G^~\ 
whereas for larger V values, one has G^ = A n < A c = 
G + G^-\ This results from the fact that, in the Hal¬ 
dane phase, the elementary excitations are all gapped, 
with a minimum occuring at ka = tt 1231 ES]- For a 
fixed value of 7/, corresponding to a fixed value of the 
ion anisotropy, D, in the corresponding spin Hamilto¬ 
nian, G (±) 

increases with increasing V (i.e. A), whereas 
G decreases. This can be understood by starting at the 
Heisenberg point, (D = 0, A = J), where, due to SU( 2) 
invariance of the spin Hamiltonian, all single particle ex¬ 
citation energies are the same, so that = G^\ In¬ 
creasing the ion anisotropy, D, i.e. going toward the Mott 
Phase (or decreasing A) gives rise to a smaller in-plane 
gap (i.e. the elementary charge gap) G^ < G^. 


This evolution of the neutral and charge gaps can 
be seen in the behavior of the structure factor S(k,uj). 
We emphasize that even though the structure factor 
S(k,uj) is a neutral excitation, i.e. conserves the total 
number of bosons, it also couples to the two-particle 
continuum composed of elementary charge excitations 
e^\k — q) + e^~\q). As explained above, in the limit 
k 0, the minimum energy corresponds to the charge 
gap A c = G( + )+G( - ), such that even if we expect S(fc, uS) 
to vanish at ka = 0, the value for A c can be obtained by 
extrapolating the behavior of S(k,uj) around ka = 0. 


In the Mott phase, top Fig. [2j one can clearly see that 
the gap at ka = tt is much larger than the gap at ka = 0. 
Note that the gap ka = tt differs from G^\ since, in 
the MI phase, e^\k) is minimum at ka = 0 and maxi¬ 
mum at ka = tt. The neutral and the charge gap have 
the same value G^ + ^ + G^~\ which can be obtained from 
S(k, id) by extrapolating the gap value to ka = 0. At the 
Mott-Haldane transition, bottom Fig. [2j S(k,ur) exhibits 
(almost) gapless excitation around ka = 0, whereas the 
excitation is cleary gapped at ka = tt. Since the transi¬ 
tion corresponds to breaking both hidden Z 2 symmetries, 
both the neutral and the charge gaps vanish, correspond¬ 
ing to vanishing elementary charge excitations gap , 
but a finite elementary neutral excitation gap G^. 
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Figure 2. TEBD excitation spectra. TOP: V/U = 0.4, 
t/U = 0.25 MI phase. The gap at ka = tv is much larger 
than the gap at ka = 0. A c = A n — S(k —* 0,u;). The white 
dashed line in the right panel shows the value of the gap (small 
but non-vanishing). BOTTOM: V/U = 0.58, t/U = 0.25 at 
the MI-HI transition. The excitations are almost gapless near 
ka — 0, whereas the excitation is cleary gapped at ka = n. 
See text. 



Figure 3. TEBD excitation spectra inside the HI phase. 
TOP: V/U — 0.7, t/U — 0.25, the gap at ka = i r, i.e. 
G^\ decreased while the gap at ka — 0 is nonvanishing, 
A n = A c = + G^ / 0. The white dashed line gives 

the value of the gap. MIDDLE: V/U = 0.75, t/U = 0.25 
inside the HI phase at the symmetric point. The gap at 
ka = 7r, G^°\ has the same value as the gap at ka = 0, 
G^ + G^~\ indicated by the white dashed line. BOTTOM: 
V/U — 0.79, t/U — 0.25 inside the HI phase. The gap at 
ka — 7T, A n = G (0) , is smaller than the charge gap at ka = 0, 
A c = + G^~\ The white dashed line in the left plot cor¬ 

responds A n = G^°\ whereas, in the right plot, it corresponds 
to A c = G m +G { -\ 


Inside the Haldane phase, but for a value V < V c , 
top Fig. [ 3 J we see that the gap at ka = 7r, i.e. G^\ has 
decreased whereas the gap at ka = 0, i.e. A n = A c = 
GW + is non vanishing. At the value V ~ V c , 

middle Fig. [ 3 J the gap at ka = 7 r, G^\ has almost the 
same value as the gap at ka = 0, G^ + G^~\ When 
V > V c , bottom Fig. [3j we see that one is now in the 
opposite situation: A n = G^ is clearly smaller than 
A c = gW+GH, 

At the transition between the Haldane phase and the 
charge density wave phase, top Fig. [ 4 J A n = G^ van¬ 
ishes while the charge gap at ka = 0, G^ + G^ re¬ 
mains finite. This corresponds to the fact that across 
the transition, the hidden Z 2 symmetry along the z-axis 
remains broken, whereas the hidden Z 2 symmetry along 
in the XY plane is restored. Since the string order O z 
is broken by charge excitations but left invariant under 
neutral excitations, i.e. commutes with the Sf operators, 
only the charge gap is protected by the finite value of the 
order O z and remains finite at the transition. On the 
other hand, since the string order O x , which is broken 
by neutral excitations (i.e. applying ), vanishes at the 
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transition, the neutral gap has to close at the transition. 


Finally, in the CDW phase, bottom Fig. [4] both the 
neutral and the charge gap increase, but still having 
A n = G<°> < GO + GO = A c . 




Figure 4. TEBD excitation spectra. TOP: V/U = 0.82, 
t/U = 0.25 at the HI-CDW transition, At ka — i r, A n = 
— 0 whereas A c — G^ + G^ at ka =■ 0, remains fi¬ 
nite and is given by the white dashed line in the right plot. 
BOTTOM: V/U = 0.84, t/U = 0.25 in the CDW phase, the 
neutral gap at ka = n is finite, but has smaller value than the 
charge gap: A n = G ^ < G^ + G ^ = A c . 
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V 




Figure 5. Excitation spectra obtained when adding a boson 
to the system, i.e. corresponding to the operators Ai — bi 
and Bj = in Eq. §. In both plots, the vertical offset 
corresponds to the chemical potential for adding a boson. 
The parameters for the top plot are the same as in Fig. [ 5 ] (top 
panels), i.e. the Mott Phase. One clearly see that the mini¬ 
mum of the charge excitation is obtained at ka = 0, the value 
at ka = 7r, being much larger. The bottom plot corresponds 
to bottom panels Fig. [3j i.e. the Haldane phase where the 
neutral and the charge gap are different. The lowest excita¬ 
tion occurs at ka = ty. Around ka — 0, one has a two-particle 
continuum, made of one neutral excitation and one charge 
excitation. 


Finally, in Fig. [5j we show the excitation spectrum 
when adding a boson to the system, i.e. corresponding 
to the operators Ai = bi and Bj = fet in Eq. (§. In both 
plots, the large vertical offset corresponds to the chemical 
potential for adding a boson, /i + . By definition, the value 
of the charge gap is — fi ~, where is the chemical 
potential for removing a boson. Therefore the minimum 
of the excitation spectrum can be written £l+ A c /2 , where 
fl = (/i + + fi ~)/2 is the average chemical potential. The 
parameters for the top plot are the same as in Fig. [2] (top 
panels), i.e. the Mott Phase. One clearly sees that the 
minimum of the charge excitation is obtained at ka = 0, 
the value at ka = tt being much larger. The bottom plot 
corresponds to bottom panels Fig. [3j i.e. the Haldane 
phase where the neutral and charge gaps are different. 
The minimum of the excitation occurs at ka = 7r, whereas 
around ka = 0, one has a two-particle continuum, made 
of one neutral excitation and one charge excitation, the 
minimum value, at ka = 0, being fi + + G^\ i.e. A c /2 + 

GW. 


IV. SUPERSOLID PHASE 


The hallmark of the supersolid phase is the presence 
of both a long range diagonal (density) order and super¬ 
fluidity. A typical density profile, obtained using DMRG 
for U = 1, V as 0.75 and t = 0.2, is shown in Fig. [6j The 
oscillations of the density between 0.25 and 2.25, around 
the average value n = 1.25 signal long range density order 
but which, nonetheless, is not in the CDW phase since 
the average density, n = 1.25, is not commensurate. 
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Figure 6. Density profile in the supersolid phase U = 1, 
V = 0.75 and t = 0.2, n — 1.25. The density oscillates 
between, 0.25 and 2.25, indicating long range density order 
at the incommensurate average density n = 1.25. 
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Figure 8. t = 0.2, n as 1.167. As in Fig. [7| the periodicity of 
the lower part of the excitation spectrum can be understood 
from the low energy effective hamiltonian, which is an AF spin 
1/2 chain in a finite magnetic field. The additional gapless 
mode at ka = tt/3 corresponds to the low energy incommen¬ 
surate modes of the spin chain at a finite magnetization. 


Figures [7| [8} [9| and [To| show the the dispersion for sev¬ 
eral values of the doping: n = 1.25, n = 1.167, n = 1.125 
and n = 1.08333 respectively. All data were obtained for 
L = 96 sites, V = 1, V = 0.75, t = 0.2 and the system is 
in the supersolid phase. 

As expected in the supersolid phase, the system ex¬ 
hibits gapless excitations at k = 0 and k = tt, but one 
clearly sees additional gapless excitations at a momen¬ 
tum kss that depends on the density. It turns out that 
the value of kss is in excellent agreement with the value 
27R s )n (see below), where S n = n — 1, i.e. ka = n/2 
(ka = 7t/3, ka = 7t/4, ka = 7t/6 ) for 5 n = 1/4 (S n = 1/6, 
£n**l/Mn=*l/12). 



Figure 9. t = 0.2, n — 1.125. As in Fig. [8] The additional 
gapless mode at ka = 7r/4 corresponds to the low energy 
incommensurate modes of a spin-1/2 chain at a finite magne¬ 
tization. 



Figure 7. t = 0.2, n=1.25. The gapless excitations at ka = 0 
emphasize the superfluid nature of the phase. The lower part 
of the excitation spectrum has a periodicity ir/a reflecting the 
2 a periodicity of the low energy effective hamiltonian, which is 
an AF spin 1/2 chain (see text). The additional gapless mode 
at ka = 0 . 57 T corresponds to the gapless mode at 2 ka = tv of 
the effective AF chain. 



Figure 10. t — 0.2, n = 1.08333. As in Fig. [8j The ad¬ 
ditional gapless mode at ka = ty/6 corresponds to the low 
energy incommensurate modes of a spin-1/2 chain at a finite 
magnetization. 
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A. Effective Spin-1/2 Heisenberg model of the 
supersolid 


We consider the situation where the supersolid phase 
has a density n = 5/4, see Fig. [6j This occurs when one 
dopes the system in the CDW phase, i.e. for V > U/ 2. 
The density pattern, in the limit U £, obtained from 
both the DMRG and the QMC computations show that 
the ground state has a (nearly) vanishing density on al¬ 
ternate sites and that the other states are (almost) built 
on either the n = 2 or n = 3 Fock states. We, therefore, 
expect the low energy excitation to be given by an effec¬ 
tive spin-1/2 Heisenberg model, where we map |3) (|2)) 
to | t) (| I)), and remove the state with vanishing density. 
The new chain has therefore an effective lattice spacing 
equal to 2 a. The effective interaction arises from the 
virtual hopping of the bosons to the empty sites. Taking 
into account the different intermediate states, one obtains 
the following effective spin-half Heisenberg hamiltonian: 


Heft = Y, ^(StSr +1 +SrS-+i +1 )+\ e gS?S? +1 -B eB S?, 

i 

(8) 

where 

Jeff = (41/- 2 u ) 

A 0,2 ( 3 2 3 _ 2 _\ 

eff \W-2U AV-U 5V-2 U 3 V-UJ 

B e ft = 2,5 p + 2 t ( 5F _ 2C / “ 3 V-u)' 

( 9 ) 

where Sp = p — 2U. A e ff is positive for a large range of 
(V,U) values, so that the preceding Hamiltonian corre¬ 
sponds essentially the AF spin-half Heisenberg model in 
a magnetic field (the negative sign of J e ff can be removed 
through the mapping ( S x ,S y ) -+ — (S x ,S y )). 

Defining the ratio p = 2 V/U, one obtains: 


Jeff = J u Gyh) 

_ 2t 2 / 3 2 6 4 \ 

eff “ U \2p-2 + 2p- 1 5,9-4 3p-2j 

B -« = 2 ^+ 4(2 (+-+)• 

( 10 ) 

For B e ff = 0, the system is in the AF ( XY ) phase when 
A e ff < | Jeff |, (A e ff > | J e ff |) [MICH]. From the preced¬ 
ing expressions, the ratio A = A e ff/| J e ff| starts from the 
value A = 2 at p = 1 and then decreases. The isotropic 
point A = 1 is crossed around p « 1.15, such that for 
V = 0.75 U, i.e. p = 1.5, the system is the XY phase. 
The ground state has therefore a vanishing magnetiza¬ 
tion, corresponding to an average density n = 5/4. For 
t = 0, this corresponds to Sp = 0, i.e. to the boundary 
between the n — 1 and n — 3/2 CDW. Then, for any fi¬ 
nite B e ff, the average magnetization is positive (negative) 


corresponding to an average density larger (less) than 
5/4. In addition, the effective spin correlations exhibit 
spatial oscillations whose period depends on the magne¬ 
tization, which, in turn, leads to gapless excitations at 
finite momentum. More precisely, starting from the cor¬ 
relation functions of the spin-1/2 chain obtained using 
the bosonisation approach [37] and taking into account 
that J e ff is negative and that the lattice spacing is 2a, one 
can show that both the in-plane ((S+S~)) and the out-of¬ 
plane (( S Z S Z )) correlations yield oscillations correspond¬ 
ing to a wavevector k^a = it( m + 1/2), where m is the 
magnetization. The latter is related to the average den¬ 
sity n of the Bose-Hubbard model: m = 2n — 5/2, such 
that rn + 1/2 = 2(n — 1) = 2 5n, where Sn = n — 1. There¬ 
fore, the gapless excitations correspond to koa = 2ir5n, 
in perfect agreement with Figs. [7J [8j [9] and [To] 

Finally, one can see in Figs. [7| [81 [9] and 10 that the 


effective period of the lowest part of the excitation spec¬ 
trum is 7r/2a, corresponding the doubling of the lattice 
spacing. This is not true for the higher excitations which 
are gapped and most likely involving the empty sites. 


B. SS-SF transition 


We have also studied the evolution of the structure 
factor S(k,u) across the supersolid-superfluid transition, 
at fixed density and interaction strengths, increasing the 
hopping amplitude from t/U = 0.24 (SS) to t/U = 0.3 
(SF). 


1 . TEBD results 


From the density plots, see Fig. [TT] (top), Fig. [I2| (top) 
and Fig. 13 (top), we see that the SS-SF transition is 
driven by the disappearance of the spatial modulation. 
This behavior is also predicted by the standard mean- 
field theory where the ground state of the system is as¬ 
sumed to be a tensor product of onsite wavefunctions (see 
below). 


The disappearance of the spatial modulation results in 
an opening of the gap at ka = 7 r, see Fig. 13 (bottom), 
which is well described by the mean-field theory, since it 
only amounts to a change of the spatial periodicity, i.e. 
from 2 a to a, of the effective Hamiltonian [38] . 


On the other hand this simple mean-field cannot cap¬ 
ture the long range quantum correlations that lead to the 
gapless modes at ka = kss& = ^irdn and the mapping 
to the spin-1/2 is no longer valid close to the transition 
since one cannot neglect previously empty sites. From 
that point of view, the exact fate of these gapless mode 
is still lacking a physical explanation. 
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Figure 11. The SS phase at U = 1, V = 0.75 and t = 0.24, 
n — 1.25. Top: Density profile: The CDW order is still 
almost perfect. Bottom: The gapless modes at ka = 0.57T are 
still visible, but their contributions to S(k,u) have a smaller 
weight when compared with Fig. [ 7 ] 


Figure 13. The superfluid phase at U = 1, V = 0.75 and 
t = 0.3, n m 1.25. Top: The density profile no longer shows 
CDW pattern. Bottom: The only gapless mode is at ka — 0, 
as expected in the SF phase; at ka = i r, the system is now 
gapped, 




kafir 


Figure 12. The supersolid phase at 77 = 1, V = 0.75 and t = 
0.26, n = 1.25. Top: The density pattern does not show a well 
defined CDW. Bottom: The gapless mode at ka = 0.57T has 
almost disappeared, but the system is still gapless at ka = ir 


2. Mean-field results 

As explained above, a well known mean-field method to 
solve the Bose-Hubbard model is the Gutzwiller ansatz, 
where the ground state wavefunction is assumed to be a 
tensor product of onsite wavefunctions: 

N max 

W = G)l^> where | ipi) = ^ f n ,i\n,i). (11) 

i n =0 

|n, i) represents the Fock state of n atoms occupying the 
site z, n max is a cut off in the maximum number of atoms 
per site, and f n j is the probability amplitude of having 
the site i occupied by n atoms. 

Minimizing the mean-field energy (T|i7|T) over the 
allows us to determine the mean-field ground 
state properties as functions of the different parame¬ 
ters (7/, £, V, fi). For instance, the superfluid phase cor¬ 
responds to a non vanishing value of the order param¬ 
eter (\k|6|\Ii), whereas the Mott phase corresponds to 
a vanishing order parameter and the ifi are pure Fock 
states. In the CDW phase, the order parameter ( b) van¬ 
ishes; the density, (n(ka = 0)) and the staggered density, 
(n(ka = 7r)), have the same value. The supersolid phase 
corresponds to non-vanishing values for both (b(ka = 0)) 
and (b(ka = 7r)); the density still exhibits oscillations at 
ka = 7r. The superfluid phase corresponds to a homo¬ 
geneous density and only the ka = 0 order parameter 
( b{ka = 0)) has a non-vanishing value. 
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We present mean-field results for U = 1, V = 1.5 
and fi = 1.8. Note that since the chemical potential 
is fixed, the density changes as t/U is changed. Fig- 
shows the different quantities as functions of t/U. 
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For 0 < t < 0.25, the system is in the CDW insulating 
phase: The order parameter ( b ) vanishes and the density 
and staggered density have the same value; the CDW 
corresponds to a density pattern • • • 2020202020 • • •. For 
0.25 < t < 1.26 the system is in the supersolid phase: 
Both (b(k = 0)) and (b(ka = i r)) are non-vanishing. The 
density still presents oscillations at ka = 7r. For t > 1.25, 
the system is in the usual superfluid phase. 

The mean field approach also allows us to compute 
the excitation spectrum. Since in both the CDW and 
SS phases the periodicity of the ground state is 2a, 
the spectrum is defined in the reduced Brillouin zone 
[—7r/2a, 7r/2a] and has two branches. In the CDW phase, 
the elementary excitations are gapped, as expected. In 
the SS phase, the lower branch becomes gapless with a 
linear behavior around k = 0. At the SS-SF transition, 
the periodicity of the ground state goes back to a, so 
that the two elementary excitation branches merge at 
ka = 7r/2. 



t/U 

Figure 14. (Color online) Mean field phase diagram for U = 
1, V = 1.5 and n — 1.8 as a function of t. For 0 < t < 
0.25, the system is in a CDW insulating phase: the order 
parameter (b) vanishes and the density and staggered density 
have the same value. The CDW corresponds to a density 
pattern • • • 2020202020 • • •. For 0.25 < t < 1.26 the system 
is in the supersolid phase: both (b(k = 0)) and (b(ka = tv)) 
are non-vanishing. The density still presents oscillations at 
ka = 7r. For t > 1.25, the system is superfluid. 


V. UNDERDOPED HALF-FILLING CDW 

In this section, we compare the structure factor ob¬ 
tained in the SS phase with the one for the phase be¬ 
tween the half-filled CDW and the superfluid phase, see 
Fig. [l5j A typical density profile is shown in Fig. [16} the 
parameters are U = 1, V = 0.75 and t = 0.1, correspond¬ 
ing to an average density n = 0.4375. One clearly sees 


that the density pattern is different from the one in the 
supersolid phase: the long wavelength modulation of the 
CDW is a signature of a vanishing DLRO, in contrast 
to the SS phase. On the other hand, one has an over¬ 
all power law decay of the ODLRO g{pc) oc x~ x ' 2K ^ but 
with a coefficient K < 1/2 emphasizing that the SF can 
be localized with a single impurity [24] . The difference 
with the SS phase also appears in the structure factor, 
Fig. [17| One has only two gapless modes, one at ka = 0 
and one at ka ~ 0.87T, but the excitations at k = tt are 
gapped. Therefore, the periodicity of the lower part of 
the spectrum is just 2ir/a and not 7r/a as in the SS phase. 


Since, at very low values of t, the density pattern for 
the half-filled CDW phase is • • • 01010101 • • •, there is a 
natural mapping onto a spin-half AF Heisenberg model 
with a vanishing total magnetization: |0) —* | and 
|1) —» | t)- The underdoped CDW phase corresponds 
then to a non vanishing total magnetization S z ~ n— 1/2. 
However, contrary to the SS phase, there is no sim¬ 
ple way to get the effective J and A coefficients: the 
initial state • • • 01000101 • • • and a state after one hop¬ 
ping ••• 01001001 •• • are actually degenerate in the limit 
t 0, thereby preventing a proper expansion of the 
Bose-Hubbard Hamiltonian. Nevertheless, one can still 
argue that the gapless mode at ka ~ 0.87T corresponds to 
the incommensurate gapless mode appearing in the spin- 
spin correlation function for the spin-half AF Heisenberg 
model in a finite magnetic field. In particular the value 
k^a = 0.87T is compatible with the bosonisation predic¬ 
tion 27m m 



Figure 15. (color online) Detail of the n — 1/2 lobe (from 
QMC) where we also determined the constant K lines for 

K= 1, 1/2. 
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Figure 16. Density profile in the underdoped 1/2 CDW phase 
U — 1, V = 0.75 and t = 0.1, n = 0.4375. The overall 
modulation of the density emphasizes the difference with the 
SS phase. 



Figure 17. (color online) Structure factor in the underdoped 
CDW: U = 1, V = 0.75 t = 0.2, the density is n = 0.40625. 
The gapless mode at ka = 0 indicates the ODLRO. The gap¬ 
less mode at ka ~ 0.87T is compatible with the bosonisation 
prediction 2tt n for an AF spin-1/2 chain in a finite magnetic 
field. 


VI. CONCLUSIONS 


In summary, we have studied the excitation spectra of 
the extended Bose-Hubbard model. Along, the MI-HI- 
CDW transition, the dynamical structure factor exhibits 
behavior similar to the spin-spin correlation for the 5=1 
Heisenberg model. For instance, it shows a difference 
between the neutral and charge gaps in the HI phase. 

In the SS phase, we have shown that the system has 
additional gapless modes at a k value that depends on 
the average density of the system. They can be mapped 
to the incommensurate gapless modes of an AF spin-1/2 
chain at finite magnetization. They are a signature of 
the DLRO present in the SS phase. These modes fade 
away when moving towards the SF phase, and, in addi¬ 
tion, a gapped mode appears at ka = 7r, marking the 
disappearance of the DLRO. 

Finally, we have shown that underdoping the CDW at 
half-filling, the excitation spectrum differs from the one 
in the SS phase, emphasizing that even though the sys¬ 
tem exhibits superfluidity and oscillations in the density, 
there is no DLRO. 
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